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ABSTRACT 


A  compilation  of  the  analytical  operations  for  the  structural  anal¬ 
ysis  of  a  truncated  spherical  shell  based  jointly  upon  an  evaluation  of 
research  effort  and  operational  performance  is  presented.  The  four 
primary  problem  areas  of  loading,  shell  conversion,  stability,  and 
strength  which  confront  the  engineer  in  the  structural  design  of  such  a 
shell  are  discussed  in  detail.  In  each  of  these  problem  areas,  com¬ 
putational  methods  or  approaches  are  offered  which  may  be  adopted  as 
design  procedure. 
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GUIDE  TO  THE  STRUCTURAL  DESIGN  OF  TRUNCATED 
SPHERICAL  RIGID  GROUND  RADOMES* 

INTRODUCTION 

The  Impetus  for  this  report  was  generated  by  the  Air  Force's  program 
in  the  development  of  truncated  spherical  rigid  radomes  for  the  environ¬ 
mental  protection  of  ground  radar  ana  communication  antennas.  As  such, 
this  report  should  be  of  particular  value  to  the  engineer  concerned  with 
the  structural  analysis  or  design  of  these  radomes;  but  the  principles 
involved  are  equally  applicable  to  any  shell  structure  of  the  configura¬ 
tion  mentioned.  Although  this  report  is  primarily  directed  toward  radome 
designers,  no  attempt  will  be  made  to  discuss  the  electronic  interaction 
of1  antenna  and  radome.  Rather,  it  will  be  assumed  that  the  radome  or 
spherical  shell  design  lias  progressed  to  the  point  where  a  tentative  type 
of  structure  has  been  selected  based  upon  transmission  characteristics; 
and  one  must  now  determine  by  analysis  what  shell  thickness,  member  sizes, 
skin  thicknesses,  and  so  forth,  are  required  for  structural  adequacy  un¬ 
der  the  imposed  loading. 

A  significant  portion  of  the  material  contained  in  this  report  was 
developed  for  Rome  Air  Development  Center  by  Armour  Research  Foundation 
under  contract  AF  30  (6i>2)  -  1R60.  Appendixes  A,  C,  and  D  of  the  final  re¬ 
port1  on  this  contract  have  been  Included  in  total  in  this  report.  This 
study  program  reflects  the  Air  Force's  belief  that  additional  theoretical 
and  experimental  work  was  ami  still  Is  required  in  this  area.  This  belief 
was  established  from  the  kaleidoscope  of  analyses  which  have  been  sub¬ 
mitted  in  response  to  operational  requirements  for  radomes.  Tills  report 
represents  an  effort  to  establish  common  ground  for  contractor  and 
Government  engineers  alike  involved  In  the  design  of  radomes. 

In  its  capacity  as  responsible  Center  for  the  development  of  ground 
radomes,  RADC  lias  had  the  opportunity  to  gather  data  on  all  phases  of  the 
program  from  basic  research  to  operational  suitability.  The  procedures 
recommended  herein  are  based  jointly  upon  the  evaluation  of  laboratory 
test  data  and  appraisals  of  rigid  radomes. 

LOADI NG 

Cased  upon  observed  field  tests  and  reports  from  operational  sites, 

it  has  been  established  that  wind  loading  is  of  primary  interei  ^  the 
♦Released  14  February  1961, 
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structural  engineer.  Although  snow  and  ice  loadings  represent  a  transmis¬ 
sion  hazard  when  one  is  considering  a  radome,  experience  has  taught  that 
for  the  typical  geome try -under.  consideration  (Figure  1)  there  is  not  suf- 
ficient  build-up  of  ice  or  snow  to  warrant  serious  structural  considera¬ 
tion.  Therefore  this  section  of  the  report  will  discuss  only  the  wind 
loading.  Its  analytical  representation  is  presented  in  Appendix  A. 

Two  types  of  analytical  representations  of  wind  load  distribution  on 
the  spherical  portion  of  the  radome  used  in  the  past  are 

.  potential  flow  distribution  on  a  complete  sphere,  and 
.  distributions  which  approximate  wind  tunnel  data  on  radomes  of 
particular  geometry. 

The  experimentally  determined  pressures  agree  well  with  potential 
theory  on  the  forward  half  of  the  dome,  but  are  substantially  less  than 
those  predicted  on  the  rear  half.  The  differences  are  due  principally  to 
the  fact  that 

.  air  does  not  act  as  an  ideal  fluid  at  the  Reynolds  numbers  in¬ 
volved, 

.  turbulent  flow  occurs  over  the  rear  half  of  the  dome, 

.  the  presence  of  the  base  geometry  (e.g.,  tower)  cannot  be  ignored. 


Figure  I.  Typical  Radome  Geometry 
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The  analytical  representations,  approximating  the  experimental  data, 
imply  a  pressure  attribution  which  is  circularly  symmetric  with  respect 
to  the  axis  containing  the  stagnation  point, and  in  the  ..direction  of  flow.  - 
The  experimental  results  of  Reference  2  on  a  3/4  tower-mounted  sphere 
(Figure  2,  Appendix  A)  show  this  assumption  of  symmetry  to  be  valid  rea¬ 
sonably  close  to  the  stagnation  point  and  also  in  the  rear  half  of  the 
dome,  but  in  error  on  intermediate  sections.  However,  a  more  complete 
analytical  representation  of  the  experimental  data  is  impractical,  es¬ 
pecially  from  the  viewpoint  of  subsequent  analytical  solution  of  the 
shell  equations. 

It  is  believed  that,  as  a  basis  for  stress  analysis,  a  general  poly¬ 
nomial  representation  of  loading,  providing  that  sufficient  terms  are  re¬ 
tained,  is  preferable  to  the  loading  as  predicted  from  potential  theory 
on  a  complete  sphere.  The  question,  then,  concerns  the  adequacy  to  which 
the  particular  series  considered  represents  the  available  experimental 
data,  say,  In  the  vertical  plane  of  symmetry  of  the  dome.  A  three-term 
series  preserves  stagnation  pressure  and  total  lift  and  drag.  As  such, 
it  closely  approximates  the  experimental  distribution  near  the  stagnation 
point,  but  deviates  markedly  in  the  rear  of  the  dome.  The  four- term 
series  evaluated  in  Appendix  A  preserves,  additionally,  the  location  of 
the  maximum  suction  pressure  (see  Figure  4,  Appendix  A),  of  course,  any 
uesirod  degree  of  accuracy  can  be  achieved  by  retaining  sufficient  terms 
in  the  series.  Subsequent  solution  of  the  shell  equations  becomes  more 
cumbersome,  but  remains  feasible.  Additional  study  is  required  to  evalu¬ 
ate  the  practical  significance  of  Including  terms  beyond  the  fourth. 

The  situation  is  less  involved  with  respect  to  stability  considera¬ 
tions.  Here  it  seems  reasonable  to  base  the  dome  design  on  the  stagnation 
pressure  assumed  to  act  uniformly  on  the  buckled  area-  Even  if  subsequent 
study  indicates  the  nonuniformity  of  the  pressure  to  be  significant,  the 
distribution  according  to  potential  theory,  being  valid  for  almost  the 
entire  forward  half  of  the  dome,  should  be  adequate. 

SHELL  CONVERSION 

In  the  theoretical  work  discussed  in  the  following  sections  of  this  re¬ 
port,  all  formulas  and  computations  will  be  based  upon  a  structure  with  a 
uniform  wall  thickness  of  homogeneous ,  isotropic  material.  Unfortunately, 
the  engineer  rarely  finds  himself  confronted  with  such  an  idealized 
structure.  More  often,  the  shell  is  a  space  frame  with  orthotropically 
oriented  members  or  of  sandwich  construction  with  a  wall  of  high  density 
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high  strength  skins  and  low  density,  low  strength  core.  To  reduce  these 
types  of  shells  to  a  form  which  can  be  handled  by  a  general  analysis,  they 
are  converted  to  an  equivalent  homogeneous,  isotropic,  uniform  thickness 
shell  of  suoh  modulus  of  elasticity  and  thickness  as  to  preserve  the  ex- 
tensional  and  bending  stiffness  of  the  original  structure.  The  elastic- 
ally  similar  shell  thus  created  will  develop  a  loading  and  deflection  pat¬ 
tern  substantially  identical  to  that  of  the  original.  The  moments, 
shears,  and  thrusts  computed  on  the  basis  Of  the  equivalent  shell  are 
then  considered  as  applied  to  the  actual  shell  and  existing  stress  levels 
computed. 

Appendix  B  presents  methods  for  the  conversion  of  a  sandwich  structure 
and  a  space  frame  structure  to  their  elastically  equivalent  uniform  shells, 
GENERAL  STABILITY 

Depending  on  construction,  primary  instability  of  a  rigid  radome  shell 
can  be  caused 

.  by  "oil  canning"  of  the  stagnation  cap, 

.  by  local  buckling  of  the  stiffeners  in  the  case  of  space  frames, 
or 

.  by  local  buckling  of  facings  on  one  or  more  panels  of  sandwich 
construction. 

To  establish  the  general  stability  of  the  structure  of  interest,  tlm 
designer  must  first  convert  the  structure  to  its  homogeneous ,  uniformly 
thick,  elastically  equivalent  shell.  This  procedure  is  outlined  in  Appen¬ 
dix  II.  With  the  structural  characteristics  of  tl.e  equivalent  shell  thus 
established,  one  may  proceed  to  determine  the  critical  buckling  pressure 
or  general  stability  of  the  structure. 

As  a  part  of  this  report,  the  general  stability  of  a  homogeneous 
spherical  shell  has  been  considered  on  the  basis  of  energy  techniques 
(Appendix  C)  from  which  the  following  has  been  indicated: 

1.  For  design,  it  is  reasonable  to  consider  the  stagnation  pres¬ 
sure  as  being  uniformly  distributed  over  the  buckled  area. 

2.  For  shells  of  the  geometry  considered  (Figure  1),  equation 
(37) ,  Appendix  C,  predicts  that  the  dimpled  area  will  intersect  the  base 
circle  for  P4.  This  means  that,  for  "thicker"  shells,  the  base  con¬ 
nections  provide  a  stiffening  efl'ect  not  considered  in  theory. 

3.  The  maximum  stagnation  or  buckling  pressure  is  given  by  equa¬ 
tion  (3fi)  ,  Appendix  C,  ana  tl.e  estimateu  size  of  dimple  by  equation  (37), 
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Appendix  C. 


uniform  pressure  is  still  somewhat  of  a  mystery..  The  problem  is  aom- 
pounded  for  a  truncated  sphere  under  nonuniform  and  especially  nonaxi sym¬ 
metric  conditions,  as  is  the  case  for  radomes.  Though  lacking  a  fundamental 
understanding  of  the  phenomena,  it  is  believed  that  the  present  findings 
contribute  to  a  rational  basis  for  design.  Continuing  study  should  lead 
to  further  improvements  of  design  procedures. 

LOCAL  STA3ILITY 


Investigation  of  the  local  stability  of  an  individual  stiffener  in  a 
space  frame  proceeds  along  the  lines  of  u  conventional  beam-column  analy¬ 
sis.  The  axial  load  per  stiffener  is  established  from  the  normal  loads 
per  unit  length  of  shell  determined  from  tho  membrane  analysis  (Appen¬ 
dix  !))  multiplied  by  the  effective  spacing  (Appendix  11)  of  the  stiffeners. 
This  uxiul  load,  along  with  the  transverse  load  transferred  from  the  mem¬ 
brane,  forms  the  combined  loading  picture  on  the  stiffener.  The  aegree 
of  support  or  "end  fixity"  assigned  to  tho  stiffener  is  a  direct  function 
of  the  construction  details  peculiar  Lo  the  structure  of  Interest.  Its 
evaluation  depends  upon  the  designer's  experience  and  a  thorough  knowledge 
of  these  details,  and  uocs  not  lonu  itself  to  general  treatment.  With  the 
loading  and  boundary  conditions  established,  one  may  now  proceed  with  any 
one  of  a  number  of  familiar  methods  for  the  solution  of  indeterminate 
structures. 


In  the  case  of  sandwich  structures,  once  the  general  stability  has 
been  established,  the  local  buckling  of  the  faces  on  one  or  more  panels 
may  be  determined  by  methods  presented  in  Reference  3  of  this  report. 
STRENGTHS 


In  general,  stresses  In  radomo  shells  have  been  determined  employing 
curved  column,  simple  beam,  or  membrane  analyses,  or  combinations  thereof. 
The  curved  column  anu  beam  analyses  Ignore  the  essential  characteristics 
ol  a  shell  in  that  they  disregard  the  influence  of  continuity  ana  curva¬ 
ture  on  the  geometrical  relations  between  displacements  and  strains,  and 
on  the  equilibrium  conditions.  Membrane  theory ,  while  Incorporating  these 
effects,  takes  no  account  of  the  flexural  resistance  of  the  shell  and,  as 
such,  has  definite  limitations  In  Its  application.  As  will  be  discussed, 
these  limitations,  at  the  very  least,  make  questionable  estimates  of  the 
conditions  of  stress  in  the  vicinity  of  the  radonie  edge  constraints. 
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What  ia  required,  then,  is  the  application  of  a  general  tending- theory 
to  the  probleifl  of  a  truncated  spherical  shell.  As  of  this  writing,  suoh 
a  theory  has  not  been  developed;  however,  it  is  nonetheless  feasible,  A 
membrane  solution,  as  discussed  in  Appendix  D,  may  be  viewed  as  a  first 
step  in  obtaining  the  bending  solution. 

Membrane  theory  results  from  a  simplification  of  the  usual  small  de¬ 
flection  bending  theory  of  elastic  shells  in  which  the  moment  resistance 
of  the  shell  cross  section  is  neglected.  As  a  result,  the  state  of  stress 
at  a  point  in  the  shell  is  completely  specified  in  terras  of  two  components 
of  the  resultant  normal  stress  on  the  cross  section  and  the  resultant 
shear.  Moreover,  these  stress  resultants  constitute  a  statically  deter¬ 
minant  system  of  forces,  meaning  that  the  governing  (partial  differential) 
equations  can  be  solved,  in  principle,  without  regard  to  shell  deforma¬ 
tions.  This  has  an  important  corollary,  namely,  that  a  membrane  stress 
solution,  in  general,  cannot  be  found  which  satisfies  an  arbitrarily  spec¬ 
ified  system  of  boundary  tractions  or  displacements. 

With  application  to  rigid  radomes,  this  means  that  a  membrane  solution 
will  not  properly  describe  the  built-in  conditions  at  the  base  circle. 
Rather,  the  theory  will  predict  a  nonzero  system  of  base  displacements  and 
associated  curvatures.  Since  the  actual  edge  conditions  may  be  thought  to 
result  from  the  superposition  of  suitable  euge  shears  and  mocuenLs,  it  is 
clear  that  the  membrane  solution  will  be  essentially  meaningless  insofar 
as  base  stresses  are  concerned. 

One  approach  to  the  bending  theory  of  shells  involves  a  two-step  solu¬ 
tion.  first,  the  membrane  stresses  are  determined  for  the  fully  loaded 
shell  and  the  deformations  are  computed  at  points  of  support.  The  second 
step  involves  evaluation  of  the  bending  theory  for  the  shell  subjected 
only  to  the  necessary  moments  and  shears  at  the  points  of  support  to  re¬ 
move  the  excess  deformations  predicted  by  'he  membram  solution.  The 
complete  solution  is  then  the  superposition  of  the  stresses  determined  in 
the  two  preceding  problems.  Progress  to  date  does  not  permit  determina¬ 
tion  of  the  bending  stresses.  However,  the  membrane  solution  described  in 
Appendix  D  was  organized  so  as  to  readily  permit  the  second  step  in  the 
above  procedure  to  be  carried  out  at  a  later  date.  As  indicated  above, 
the  results  of  the  membrane  solution  are  considered  sufficiently  accurate 
for  design  purposes  at  all  points  throughout  the  shell,  except  in  the  im¬ 
mediate  vicinity  of  the  base  circle  (that  is,  the  boundary  layer  effect 
as  Riessner  terms  it  in  Reference  4) . 
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Sinoe  the  designer  ware  often  finds  himself  involved  wi±b.  other  l4>aa.' 
a.  uniform  thickness  shell  of  hOfllQgeae.QUS..m&lerlali  ihfi- approach  to  ft  V. 
stress  analysis  must  consider  .  . 

.  means  of  determining  an  "equivalent"  homogeneous  shell  to  rep¬ 
resent  the  actual  structure  (Appendix  B) , 

.  the  application  of  membrane  theory  (Appendix  D)  to  determine 
the  normal  and  shear  loads  in  the  equivalent  shell,  and 

.  the  application  of  these  loads  to  the  actual  shell  and  com¬ 
putation  of  the  resulting  stresses. 

SUMMARY 

In  this  report  and  the  appendixes  thereto,  we  have  attempted  to  es¬ 
tablish  a  norm  or  basis  for  a  more  rigorous  analytical  design  of  a  radome 
or  spherical  shell.  Beginning  witli  the  loading,  a  procedure  for  the  ana¬ 
lytical  representation  of  the  pressure  distribution  on  a  radome  resulting 
from  a  wind  loading  has  been  indicated  and  worked  out. 

Secondly,  methods  for  converting  space  frame  structures  and  sandwich 
structures  to  their  homogeneous,  uniform  wall  thickness,  elastic  equi¬ 
valent  counterpart  have  been  demonstrated.  Once  this  has  been  accom¬ 
plished,  the  designer  may  now  apply  the  "membrane  theory*  to  determine 
shell  loads  and  investigate  for  general  shell  stability. 

Third,  a  membrane  analysis,  v hich  can  be  readily  incorporated  into  a 
general  bending  theory,  is  presented  which  has  been  evaluated  for  a  four 
term  expansion  of  the  loading  and  is  considered  valid  in  itself  for  that 
portion  of  the  radome  shell  above  the  equatorial  plane.  The  results  of 
this  analysis  may  now  be  applied  to  the  actual  structure  under  considera¬ 
tion,  and  local  stability  and  strength  requirements  investigated. 

Fourth,  a  new  theory  of  the  buckling  of  a  complete  spherical  shell  un¬ 
der  uniform  external  pressure,  developed  by  Professor  L.H.  Donnell  of 
Armour  Research  Foundation,  is  presented  as  a  part  of  this  report.  An 
approximate  finite-displacement  solution  was  obtained  for  the  axlsymmetric 
buckling,  in  the  form  of  a  single  dimple  of  a  fairly  thin-walled  homogene¬ 
ous,  elastic  spherical  shell.  This  solution  provides  for  the  buckling 
pressure,  P cr,  and  the  serai-cone  angle  of  the  "dimple"  which  forms  the 
buckled  shape,  /3. 

RECOMMENDATIONS 

Based  upon  a  general  review  of  the  design  techniques,  the  following 
analytical  studies  are  recommended: 
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.  develop  and  evaluate  in  detail  a  general  bending  theory  for  uni¬ 
form,  homogeneous,  spherical  shells  subject  to  the  simplified  wind  loading 
of’  Appendix  A,  and*  incorporating'  the  membfaTiS  afrelysfs  of  Appendix  Dj  ’ ’ 

.  extend  the  Donnell  buckling  theory  presented  in  Appendix  C  to 
include  more  general  characterization  of  the  deflected  shape; 

.  application  of  a  more  exact  basis  for  determining  the  geome¬ 
trical  and  mechanical  properties  of  an  "equivalent"  homogeneous  shell. 

One  such  theory  due  to  E.  Riessner  (Reference  5)  would  appear  to  be  partic¬ 
ularly  suited  to  radome  structures. 
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APPENDIX  A 

-AHAL-Y-TICAL-R£PRE^ENTA T ION  OF  WIND  LOAD INS - - - — - “ 7 

Typical  data  on  steady  wind  pressures  incident  on  radome  shell  geo¬ 
metries  of  the  type  considered  (Figure  1)  was  obtained 'from  Reference  2, 
and  is  summarized  in  Figure  2.  With  a  view  toward  simplification,  in¬ 
spection  of  Figiire  2  suggests  replacing  the  given  isobars  by  ones  which 
are  symmetrical  with  respect  to  the  axis  of  flow  (these  isobars  would  ap¬ 
pear  as  vertical  lines  in  Figure  2).  With  this  assumption,  the  wind  pres¬ 
sure  is  a  function  only  of  the  distance  (measured  in  the  direction  of  flovO 


from  the  stagnation  point,  X  (see  Figure  3). 

For  purposes  of  stress  analysis,  it  is  found  convenient  to  express  the 
normal  pressure  as  a  finite  power  series  in  the  coordinate  X,  or,  equival¬ 
ently,  in  a  cosine  series  in  the  meridian  angle  9  (Figure  3).  The  series 
form  adopted  is  given  in  equation  (1) .  Symbols  are  defined  in  the  nomen¬ 
clature. 
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Figure  3 


Shell  Coordinate  System  and  Membrane  Forces 


(*.*) 


2  a.  (x/ft)B 

n»0 


Sa(l-  sin  oos  8)n 

n«0  n 


(1) 


2  Pn  ( <f> )  cos  ni), 
n=0 


where  a0  ■  Pg  (4>)  m  1* 

The  coefficients  Pn(<£)  can  be  expressed  In  terras  of  the  constants 
«n  by  expanding  equation  (1)  and  equating  coefficients  in  cos  n6.  The 
an,  in  turn,  are  related  to  those  features  of  the  actual  pressure  dis¬ 
tribution  which  are  to  be  retained.  For  N  *  3,  the  following  relations 
are  obtained: * 

P0  «  .10  +  ,\2  sin2  <ji, 

p  1  °  sin  4>  +  sin3  0, 


P5 


A  2  sin2  <p, 


)  sin3  4>, 


(2) 


.1 


where 


/l0  c  h0  +  (jj  +  n2  +  u3, 

J|  -  -  (nx  +  2a  2  +  Dn3) , 
A 2  c  1/2  (a2  +  3a3 ) , 

Ay  °  -  ri  3/i. 


(3) 


The  four  coefficients  Z0  .  .  .  are  determined  so  as  to  preserve  the 
total  lift  and  drag  on  the  radome  (corresponding  to  the  given  pressure 
distribution) ,  the  stagnation  pressure,  and  the  position  of  the  minimum 
average  pressure.  The  first  three  conditions  are  clearly  necessary,  where¬ 
as  the  fourth  condition  is  somewhat  arbitrary.  If  more  terms  in  equation 
(1)  were  retained,  additional  features  of  the  given  pressure  distribution 
could  be  preserved. 

In  terms  of  the  assumed  pressure  distribution,  the  total  lift,  L,  and 
drag,  U,  on  the  radome  shell  are  to  be  found  to  be: 


•^Seneral  relationships  are  derived  but  numerical  results  are  presented  only  for  H  =  3. 
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i  ■  “  Ir  jm-  «**  *■: 


i  ■  “2fl*  X  °  l"  Pn  (0)  cos  0  sin#  008  ng  dg  dtp  (4) 

n  0  0  Jt*V 


'f'O 

■  -  2 nR2  £  Pq  ( 4> )  sin  <p  cos  $ 


D  a  s  pu 8in  *  cos  8  d  area, 


D 

M  H 


(*9) 

2fi2  /  /  2  p  ($)  sin2  <p  cos  8  cos  h8  d8  d<p 

0  0  n*0  n 

<t'e. 


(5) 


B  vR2  £  p1  (<p)  sin2  4>  d(p. 
For  N  =  3,  using  equation  (2) 

L  vR 2 


'7 


(2  /10  +  Sz  sin2  <t>Q)  sin2  $0, 


b  vR2  , 

-  *  -  1(5.1.  +  12 /L)  (1  -  cos  4>0')2  (2  +  cos  <p0)  - 

q  1 5  ? 


(a) 


(7) 


-  9  lj  sin4  <}> o  oos  0O] . 

The  condition  for  the  pressure  at  the  stagnation  point,  <t>  *  rt/ 2,  0  *  0,  is 

PW  2,0)  *  , 

f7  n-0  " 

For  N  »  3  this  leads  to 

,l0  +  ;lj  +  2A?  +  4/1  ?  =  1.  (8) 

If  the  position  of  the  minimum  pressure  is  prescribed,  say,  at 
x  ■  sin  4>  cos  8  3  xn,  equation  (1)  yields 


d{P/q) 

dx 


X  =  X 


-  z  n  ( in  (1  -  Xj 


n-  1 


0. 


n=0 


(9) 


For  N  ■  3,  the  an  may  be  expressed  in  terms  of  the  An  by  means  of  equa¬ 
tion  (3) ,  and  equation  (9)  becomes: 

A1  +  4  Ai  X»  +  12  A3  =  °* 
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The  coefficients  Ax  .  .  .  A3  can  now  he  determined  from  equations  <6)  , 
(7),  (8),  and  (9)  one eL/q,  D/q,  and  X*  are  specified." 

Referenoe  8  reports  the  foUdwing_4xperime^iliy  deih rained  values 
of  total  lift  and  drag  for  a  13.8  inch  radius  radorae  at  a  stagnation 
pressure  of  q  ■  0.367  psi  (150  raph  wind,  one  atmosphere  wind  tunnel 
pressure) . 

£  -  183  lb, 

D ■»  52  lb  (extrapolated  for  full-scale  Reynolds  number). 

Using  these  values,  the  following  lift  and  drag  coefficients  were 
determined  for  a  hemispherical  radome  (</>fl  *  n/2,)  . 


®l 


L 

nE2q 


183 

(0.367)  (13.5) 2  v 


0.87, 


2  D  (2)  (52) 

TiE2q  (0.367)  (13.5)  2 tt 


Reference  2  also  indicates  that  the  minimum  average  pressure  in  the 
plane  6  «  0  occurs  at  x„  ■  sin  4>n  *  0.2  (<p%  ■  11.5°).  Using  the  above 
values  for  ,  C^,  and  X^,  the  J 's  are  found  to  be, 

A  q  *  -  1.1433 
»  -  0.0401 
42  ■  0.5466 
13  ■  0.4226. 

The  numerical  results  of  the  stress  analysis  presented  in  Appendix  D 
are  based  on  these  values.  Figures  4  and  5  show  a  comparison  of  the  aver¬ 
age  pressure  data  -'btained  from  Reference  2  and  the  above  power  series  ex¬ 
pansion  with  A  ■  3. 
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Figure  H.  Comparison  of  Pressure  Distribution 

WIND 

DIRECTION - ► 


Figure  5.  Comparison  of  Pressure  Distribution  (Polar  Plot) 


NOMENCLATURE  FOR  APFENDIX  A 

.Qn.'in  -Load  parameters. 

C0  Coefficient  of  drag, 

C£  Coefficient  of  lift, 

D  Drag, 

L  Lift, 

q  Stagnation  pressure, 

P  Pressure  normal  to  the  undeformed  shell  surface, 
pn  Pressure  distribution  parameter, 

X,  X  Shell  coordinates  (Figure  3), 

H  Radius  of  shell , 
d  Meridional  angle, 

4>  Circumferential  angle,  and 
<Pq  Angle  of  the  base  circle. 


APPENDIX  B 


SHELL  CONVERSION 


The  conversion  of  a  sandwich  or  space  frame  type  shell  to  Its  elastic 
equivalent  often  reflects  the  investigator's  Ingenuity  and  engineering 
Judgment  and,  without  a  more  exact  solution,  is  difficult  to  evaluate. 
Thus,  the  choice  of  one  method  over  another  often  is  not  unique  but  more 
a  matter  of  opinion.  The  following  methods  have  been  used  in  the  past  ami 
are  recommended  for  continued  use  until  such  time  as  further  study  yields 
a  more  rigorous  solution. 


To  create  the  elastically  equivalent  shell,  the  designer  must  preserve 
the  ex tens ional  and  bending  stiffness  of  the  actual  structure.  Represent¬ 
ing  extenslonal  stiffness  per  unit  width  of  structure  as  the  product  EA 
and  the  bending  stiffness  per  unit  width  as  the  product  El,  then  the  fol¬ 
lowing  relationships  must  be  maintained  between  the  actual  and  equivalent 
s  tructures . 


El. 


ac tu a  1 


EA 


actual 


E  1 

equ l valent 

E  ! 

'  equ l  valent 


SANDWICH  SHELL 


In  the  case  of  the  sandwich  shell,  the  development  continues  as  fol- 
iuws  (Reference  I!) : 

The  general  expression  for  the  bending  stiffness  of  a  unit  width  of 
sandwich  construction  is: 


(10) 


where 


This  expression  may  be  simplified  by  three  postulations  which  are  gen¬ 
erally  valid  in  the  case  of  radome  design: 

1.  Since  the  core  is  low  density,  low  strength  material,  its 
stiffness  may  be  neglected. 


2. 


Skin  thicknesses  are  equal  (that  is, 


tr). 
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3.  Since  the  skins  are  thin,  their  individual  stiffnesses  can  be 
neglected. _ 

Through  application  of  these  postulations,  equation  (10)  may  be  re¬ 
duced  to: 


tf{t  +  V  • 


(il) 


The  bending  stiffness  of  the  equivalent  shell  per  unit  width  is  ex¬ 
pressed  as: 


EL 


E  t  ; 
e  e 


equivalent  \ 


(12) 


Equating  the  actual  and  equivalent  bending  stiffness,  equation  (11) 
and  (12)  yield 


l'!  *5 

€  e 


3  Bf  lr 


(t  +  tc)S 


(13) 


Following  the  same  procedure  with  extensional  stiffness,  the  general 
expression  for  the  extensional  stiffness  of  a  sandwich  shell  per  unit 
widtli  is: 


EA. 


“  K,  «  -  t  )  +  t  . 


(14) 


Again  neglecting  the  modulus  of  elasticity  of  the  core  and  assuming 
equal  thickness  skins,  equation  (14)  becomes: 


^•'actual  *  A> 


(15) 


The  extensional  stiffness  per  unit  widtli  of  the  equivalent  plate  is 
expressed  as: 


^  ‘^equivalent  ’ 


(16) 


Equating  the  extensional  stiffness  per  unit  width  of  the  actual  and 
equivalent  shells,  equations  (15)  and  (1G)  yield: 


^ e  te  °  • 


(17) 
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Solving  equations  (13)  and  (17)  simultaneously  yields  the  fo Hawing 
expressions— £or--the—thlckness_-anji- modulus,  of  elasticity  of  the  eauivaient 
shell  in  terms  of  known  parameters  of  the  actual  sandwich  shell. 

/3 

t.  =  T  <«  +  *e)» 

8  .  (1.8) 

/4‘  V? 

i?  ■  —  - —  , 

*  3  t  +  tc 

With  the  thickness  and  modulus  of  elasticity  of  the  equivalent  shell, 
one  may  now  proceed  to  compute  load  distribution,  deformation,  and  buck¬ 
ling  pressure.  However,  it  must  be  kept  in  mind  that  the  actual  stresses 
in  the  dome  will  be  computed  from  the  loads,  determined  from  the  equiva¬ 
lent  shell  per  unit  width,  as  applied  to  the  actual  shell  of  interest.  In 
the  case  of  the  sandwich  shell,  the  stresses  in  the  core  and  facings  due 
to  these  loads  can  be  calculated  by  means  of  methods  given  In  Reference  3. 
SPACE  FRAME 

To  develop  the  same  equivalency  for  a  space  frame  structure  which  re- 
seniDles  an  irregular  polyhedron,  one  must  estimate  the  "effective  spacing” 
served  by  each  of  the  stiffeners  comprising  the  irregular  polyhedron.  It 
lias  been  demonstrated  (Reference  0,  p  !>5)  that  as  the  number  of  faces  of  a 
regular  polyhedron  increases,  the  stiffener  stress  under  uniform  pressure 
approaches  the  stress  in  a  sphere  of  equivalent  surface  area  similarly 
loaded  whose  thickness  has  been  determined  as  follows: 

total  vol  of  stiffeners 

t  =  — — - . _______  , 

2  x  surface  area 

Only  half  the  volume  of  the  stiffeners  has  been  equated  to  the  volume 
of  the  shell  because  of  the  orthotropic  nature  of  the  stiffener  as  opposed 
to  the  isotropic  nature  of  the  shell.  Since  the  stiffener  stress  and 
shell  stress  are  now  equal,  one  can  say  that  the  stiffener  is  now  serving 
an  effective  length  or  shell  U  s)  whose  area  is  equal  to  that  of  the  stif¬ 
fener  (  .-I  s )  , 


I _  x 

■S  S 


This  length  of  shell  will  be  used  as 
stiffeners  in  developing  the  elastically 


the  "effective  spacing"  of  the 
equivalent  shell  as  follows. 
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As  in  the  case  of  the  sandwich  shell?  the  following  equations  must  be 
satisfied!  ••• 


^actual  "  ^ equivalent  an<^ 

PA  «  &  A 

““actual  equivalent . 

The  bending  stiffness  per  unit  width  of  the  stiffened  shell  is  ex¬ 
pressed  as: 


El 


V  I 

s  s 


actual 


(19) 


With  the  bending  stiffness  of  the 
in  equation  (12)  and  equating  the  two 


equivalent  shell  again  expressed  as 
bending  stiffnesses,  we  have: 


where 


E  t  3 

e  ve 

12  K 


(so) 


Proceeding  in  the  same  fashion  witli  the  extensional  stiffness  per  unit 
wiuth: 


EA 


a  c  ( uo  1 


(SI) 


Witli  the  extensional  stiffness  of  equivalent  shell  again  expressed  as 
in  equation  (1C)  and  equuting  the  two  extensional  s tiffnesses,  we  have: 


E  A 
s  s 


E  t  , 

e  e 


(88) 


Proceeding  in  a  fashion  identical  to  that  with  a  sandwich  and  solving 
equation  (20)  and  (22)  simul taneous lj  fields  the  following  expressions  for 
the  thickness  and  modulus  of  elasticity  of  the  equivalent  shell  in  terms 
of  known  parameters  of  the  stiffened  shell, 
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Once  again  the  load  distribution,  deformations  and  buckling  pressure 
may  be  computed  based  upon  the  elastically  equivalent  shell.  To  apply 
the  load  distribution,  as  determined  from  the  membrane  analysis  (Appen¬ 
dix  H)  to  the  stiffened  shell  in  terms  of  axial  load  per  stiffener,  one 
considers  the  stiffeners  as  assuming  an  axial  load  equal  to  the  load 
carried  by  a  length  of  the  equivalent  shell  equal  to  the  "equivalent 
spacing"  developed  for  the  stiffeners.  This  axial  load,  plus  the  trans¬ 
verse  load  resulting  from  the  transfer  of  the  wind  load  to  the  stiffener 
by  the  membrane,  form  the  combined  loading  picture  under  which  the  stif¬ 
fener  is  analyzed  for  local  stability  and  strength. 


I 

I 
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NOMENCLATURE  FOR  APPENDIX  B 


^  ^acttto 1 

equi valent 
^'^actva  l 
AA e quiva lent 


b 

b 


t 

b 

in  £  n 
l  j  f. 

tf, 

Hr 

He 


t 


e 


.4 

Z 

£ 

/ 


s 

s 

a 

s 


Bending  stiffness  per  unit  length  of  actual  shell  under  con- 
sideration 

Bending  stiffness  per  unit  length  of  elastically  equivalent 
shell 

Extensional  stiffness  per  unit  length  of  actual  shell  under 
cons ideration 

Extensional  stiffness  per  unit  length  of  elastically  equiv¬ 
alent  shell 

Young's  modulus  of  elasticity  of  sandwich  skin  material 

Young's  modulus  of  elasticity  of  sandwich  core  material 

Over-all  sandwich  thickness 
Core  thickness 

Individual  skin  thickness  when  skins  are  unequal 

Individual  skin  thickness  when  skins  are  equal 
Poisson's  ratio  of  the  skin 

Poisson's  ratio  of  the  core 

Young's  modulus  of  elasticity  of  elastically  equivalent  shell 
Thickness  of  elastically  equivalent  shell 
Thickness  of  equivalently  stressed  shell 
Cross  sectional  stiffener  area 

"Effective  spacing"  of  Indivluual  stiffeners  in  a  grid  work 
Young's  modulus  of  elasticity  of  stiffener  material 
Moment  of  Inertia  of  stiffener  cross  section 
Poisson's  ratio  of  stiffener  material 

(1  -  M2) 


P 


1- 


A'c  kr 
T^kc 
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APPENDIX  C 

EHERQY  THEOB ¥_&f_JMU_BUCKUM$  (Rgfar.ftnoe  1) — 1 —  -  — . 

Introduction 

This  Appendix  presents  an  approximate  solution,  on  tshe  basis  of 
energy  techniques,  for  the  buckling  of  a  spherical  shell  under  uniform 
external  pressure.  More  specifically,  a  finite-displacement  solution  is 
obtained  for  the  ax is y mme trie  buckling,  in  the  form  of  a  single  dimple, 
of  a  fairly  thin-walled,  homogeneous,  elastic,  spherical  shell.  The  solu¬ 
tion  indicates  that  the  pressure  distribution  outside  of  the  dimple  area 
is  unimportant  insofar  as  the  buckling  pressure  is  concerned.  Hence,  it 
appears  that  for  buckling  under  wind  pressures,  the  effective  uniform  pres¬ 
sure  could  be  taken  as  the  pressure  at  the  stagnation  point  without  great 
error.  The  theory  provides  explicit  forms  for  the  buckling  pressure  and 
the  size  of  dimple  in  terms  of  the  geometric  and  mechanical  properties  of 
the  shell. 

The  present  solution  differs  from  previous  studies  of  this  type  (e.g., 
Reference  7)  in  two  principal  respects: 

.  it  takes  Into  account  a  more  general  class  of  functions  as 
characterizing  the  buckled  shape;  and 

.  it  accounts  for  initial  imperfect  lulus  in  the  geometry  of  the 

shell . 

The  theory,  of  course,  suffers  none  of  the  drawbacks  of  shallow  shell 
theory  especially  in  that  no  assumptions  as  to  the  boundary  conditions  at 
the  edge  of  the  dimple  are  Introduced. 

The  following  sections  present  the  derivation  of  the  governing  equa¬ 
tions,  the  method  of  solution,  and  the  results  obtained. 

DERIVATION  OF  EQUATIONS 

Figure  6  shows  a  section  of  a  spherical  shell  of  radius  k  and  constant 
thickness  t.  The  position  of  a  point  in  the  middle  surface  is  located  by 
means  of  the  spherical  coordinates,  r,6.  The  coordinate  Z  measures  the 
distance  of  a  point  from  the  middle  surface.  The  deformation  of  the  mid¬ 
dle  surface  is  assumed  to  be  radially  symmetric,  and  is  given  in  terms  of 
a  tangential  displacement  l>  and  radial  displacement  w,  directed  positive 
as  shown  in  Figure  6.  The  middle  surface  is  also  assumed  to  have  an  ini¬ 
tial  radial  deformation  prior  to  application  of  the  uniform  normal  pres¬ 
sure  p. 
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the  initial  radial  deviation  from  spherical  shape,  «{,  is  considered  which 
Is  the  same  shape  as  the  deflection  under  load  &>. _ .  ...  _ _  .  .. 

The  total  potential  energy  e  is  .  .. 


77  rlE 

~ 7* 


t/2 


"  ~r)iZ  f  U,!  +  £,!  *  Sve*  £*!  31,1  6,10  ‘ 


277  rz  P  co  sin  0  <20 


(26) 


For  simplicity  In  the  derivation  of  equation  (26) ,  it  has  been  assumed 
that  the  pressure  p  remains  constant  during  displacement  of  the  shell.  It 
is  believed  that  this  will  not  introduce  significant  error  if  p  is  con¬ 
sidered  to  be  the  average  pressure  during  buckling. 

The  displacements  cn  and  V  are  assumed  to  have  the  shape, 


-y 


a I  "  te  £  V  yn 


m  »  0,  2,  4  .  .  .  even  (27) 


t  2  5  ® 

p  =  —  e  ^  K- 
r 


n«l 


+  1 


It  is  seen  that  the  coefficient  fr'0  determines  tjie  amplitude  of  the  dimple, 
and  the  other  coefficients  its  precise  shape.  The  number  k  determines  the 
effective  size  of  the  dimple.  A  semi-angle  of  the  dimple,  /?,  may  be  de¬ 
fined  somewhat  arbitrarily  in  terms  of  k  for  given  r/t.  As  shown  in  Fig¬ 
ure  7,  it  seems  reasonable  to  associate  the  value  of  v  *  2.5  with  the  an¬ 
gle  /3, 

It  will  be  observed  that  the  expressions  for  displacement,  equation 
(27) ,  do  not  satisfy  exactly  continuity  conditions  on  the  opposite  side  of 
the  sphere,  6  =  180  degrees.  Tills  is  considered  to  be  of  no  significance, 
however,  since  the  displacements  are  completely  negligible  there. 

METHOD  OF  SOLUTION 

In  the  solution  presented  herein,  only  tiie  coefficients  V0 ,  Fj,  and  Fj 
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Figure  7.  Determination  of  Semi-dimple  Angle  fi 


In  equation  (27)  were  considered,  all  other  coefficients  being  taken  zepo. 
By  means  of  equations  (24)  ,  (2C>)  ,  and  (27)  ,  and  assuming  the  shell  to  be 

r  t  * 

sufficiently  thin  so  that  -  ft  »  1  and  -  ft  <<  1  ,  one  finally  obtains 

t  T ' 


3(1  -  v 2)  r 
2  TT  £  t  tt 


e 


2  +  l\  V,  + 


V)  +  u, 


G*0- 


-  »0  (V1  +  V3)  +  (2ft  ff02  +  A'0A) 


(20) 


3 

+ - 

128ft 


^The  parameter  ft  is  of  order  unity. 


1 

(2ft  2  +  A;0A)  I 


3PA'0 

2ft 


I 


—  (l~vz)Prz 

P  H  - 

.  “  -■  “  TTt*  ;  ““  '  "  ‘  '";'  ‘  ' 

A  Ur 

A  ■  .  (30) 

l 

Following  Reference  8,  the  number  U  is  considered  to  be  an  "unevenness  fao 
tor,"  defined  so  that  the  initial  radial  deviation  may  be  expressed  as 


■  ai  t  e“  y 

and  <ai> 

a  l*  Ur  A_ 

J7  "  Ut  ‘  4fr  • 


Here,  the  term  l 


is  the  approximate  half  wavelength  of  the 


radial  displacement,  taking  the  whole  dimple  as  a  full  wave.  This  expres¬ 
sion  is  similar  to  those  used  in  studying  the  effects  of  imperfections  in 
other  buckling  problems,  see  for  example  Reference  8.  In  these  previous 
applications  ,  it  was  argued  that  V  should  be  roughly  independent  of  the 
dimensions  of  the  structure  and  should  depend  only  on  the  manufacturing 
process.  The  factor  U  was  found  to  be  of  the  order  of  0.0003,  thus  a 

A  Ur 

value  of  A  =  —  ■  1  corresponds  to  an  extremely  irregular  shell  unless 
r/ t  is  very  large. 


Be  Be 

Setting  —  »  —  °  0  in  equation  (29)  according  to  the  principle  of 

virtual  work,  there  results. 


3k 

T 


(2  Pj  +  V3)  -  k0  -  0, 


3k 

~ 


4.  Of/  > 
■ 
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(2k  k'Q  Z  +  h'0  A)  a  0. 


(32) 


Solving  equation  (32)  for  and  V-}  we  find 


I 

I 

i 


r>  ■  1*  l1  *  S  <8*  *•  *  a,l 


^1,-T 

9k  [  37 


(3k  kn  +  A)  , 


Bg  3g 

Now  setting  —  •  — •  =  0  in  equation  (29) ,  and  using  equation  (33)  to 

OA  O  H'q 

eliminate  and  Fj,  there  finally  results  two  equations  relating  P,  and 
H0  that  is, 

P  *>  0.37  hi0  (1  -  0.283A  +  0.024A2  -  OASk2)  -  0,037k2  U<?,  (34) 


P  “  0.74  Vn  (1  -  0.0283A  +  0.024A2  +  0.45k2)  -  0.625 


(1  -  0.173A)  k  kQ2  +  0. 146ft2  <t'03. 

Eliminating  P  between  these  expressions,  dividing  by  h'0 ,  ana  simplifying, 
we  obtain 


(1.35  i  0.494  ^7)  k2  -  1.09  .VQ  (1-  0.173A)A>+  (1  -  0.283A  +  0.024A2)  =0.  (35) 

• 

The  desired  relationship  between  P  and  V0  (load-central  deflection  J 

curve)  can  now  be  obtained  by  assuming  values  of  A  and  solving  equation 
(33)  fork,  and  then  solving  either  of  equations  (34)  for  P.  The  results  of 
sueli  calculations  are  shown  in  Figure  8  for  the  extreme  values  of  the  uneven¬ 
ness  parameter  A  •  0  and  1.  It  is  observed  that  the  load-deflection 
curves  are  not  particularly  dependent  upon  the  value  of  A;  the  range 
0<  A  <1  being  extreme,  witli  shells  of  practical  interest  in  the  vicinity 
of  A  ■  0.  For  this  case,  the  value  P  *  0.G2  may  be  taken  as  the  nondimen- 
sional  buckling  load  (Figure  8).  From  equation  (30),  the  buckling  pres¬ 
sure  Pcr  then  has  the  value 


0.62 Kt2 

-  • 

(l  -  v2)rz 


(36) 


It  is  found  that  the  values  of  k  ao  not  vary  greatly,  and  average 
about  0.5  for  the  near  horizontal  parts  of  the  curves  in  Figure  8.  If,  as 
previously  mentioned  in  connection  witli  Figure  7  the  semi-angle  of  the 


dimple  f3  is  taken  as  7  e  2.5,  then  by  equation  (28)  witli  k  =  0.5, 
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= 3.5  j  —  (radians). 


p 


Figure  8.  Load-center  Deflection  Curves 
DISCUSSION  OF  RESULTS 

The  essential  results  of  this  stud*  are  contained  in  equations  (36) 
and  (37),  narael*,  that  the  uniform  pressure  required  to  produce  a  dimpl 
in  the  shell,  , 


is  aDout 


and  that  the  dimple  formed  has  a  semi-angle  fi  (Figure  7)  of  about 


fi 


radians 


200  /-  degrees. 

/  r 


The  theory  also  points  to  a  relatively  small  influence  of  geometric  ir¬ 
regularities  in  the  shape  of  the  shell.  Finally,  equation  (37)  shows  that 
for  reasonably  thin  shells,  the  assumptions  of  shallow  shell  theory ,  aside 
from  edge  conditions,  should  be  valid. 

The  general  form  of  equation  (36)  can  be  considered  to  be  well-estab¬ 
lished;  the  magnitude  of  the  numerical  coefficient  0.G2,  of  course,  is 
less  well-es tablished.  In  the  one  test  result  available  in  the  literature 
for  a  full  sphere,  Reference  9,  E  «  14.5  x  1U6  psi,  r/t  «  900  and  the  nor- 
pr 

mal  stress  was  —  *  2480.  These  data  correspond  to  a  numerical  coeffi- 
2t 

cient  of  0.31,  .lust  half  the  value  given  in  equation  (30).  The  material 
used  is  described  as  copper,  so  that  yielding  may  have  played  some  part  in 
its  failure.  The  value  of  fi  In  this  test  is  reported  as  eight  degrees, 
while  equation  (97)  gives  6, 7  degrees  for  this  ease,  However,  the  loca¬ 
tion  of  the  "edge"  of  such  a  dimple  Is  obviously  rather  iiiuef.1  ill  tc 

In  a  recent  series  of  unpublished  tests  at  the  Armour  Research  Founda¬ 
tion  using  ordinary  ping  pong  balls,  clear  indication  of  the  dimple  forma¬ 
tion  was  observed.  Ping  pong  balls  are  found  to  be  reasonably  spherical 
shells  for  which  r/ t  ■  50  and  E  «  5  x  105  psi  and  which  collapse  at  about 
100  psi.  Using  these  data,  the  coefficient  In  equation  (3G)  is  determined 
to  be  0.45  rather  than  0.62.  The  angle  fi  was  measured  as  30  degrees  and 
predicted  from  equation  (37)  to  be  28  degrees.  Finally,  the  experimental 
data  on  clamped  shallow  shells  reported  in  Reference  10  shows  a  coeffi¬ 
cient  of  0.4  in  equation  (36)  to  represent  a  lower  bound  on  the  collapse 
pressure;  the  value  0.62  falls  In  about  the  middle  of  the  experimental 
data. 

The  present  theory  is  approximate  on  several  counts,  chiefly  in  that 
the  buckled  shape  is  characterized  by  a  minimum  number  of  terms.  The 
inclusion  of  additional  terms  would  have  the  effect  of  reducing  the  numer¬ 
ical  coefficient  in  equation  (38).  It  is  strongly  recommended  that  work 
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(in  this  theory  be  pursued.  For  the  present,  however,  there  is  no  basis  ' 

equation  (36)  serves  as  a  realistic  basis  for  design  with  probably  a  mini¬ 
mum  of  conservation,  and  that  the  size  of  the  dimple  is  given  by  equa¬ 
tion  (37) . 


In  equations  (37)  anti  (38) ,  E  and  t  would  be  replaced  by  the  expres¬ 
sions  for  Ee  and  te,  as  presented  in  Appendix  B,  to  evaluate  the  criti¬ 
cal  buckling  pressure  for  a  space  frame  or  sandwich  design. 
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APPENDIX  D 

In  this  section*  a  stress  and  deflection  analysis  of  a  homogeneous 
spherical  radome^ shell  of  arbitrary  base  angle  and  subject  to  the  . 

pressure  loading  developed  in  Appendix  A  is  carried  out  on  the  basis  “of 
membrane  shell  theory.  It  is  intended  that  this  analysis  be  considered 
as  the  first  step  in  a  two-step  superposition  approach  to  a  general 
bending  solution.  The  second  step,  not  carried  out,  requires  application 
of  bending  theory  to  a  radome  subject  only  to  edge  tractions  sufficient 
to  compensate  for  edge  deflections  and  rotations  implied  in  the  membrane 
solution.  Toward  this  end,  the  membrane  solution  is  adjusted  to  yield 
only  radial  deflections  at  the  edge.  This  solution  is  considered  valid 
for  that  portion  of  the  radome  shell  above  the  equatorial  plane.  These 
results,  especially  with  regard  to  deflections,  are  believed  to  be  more 
complete  than  heretofore  available. 

MEMBRANE  ANALYSIS 

Stresses  -  Figure  3  shows  the . membrane  forces  which  act  on  an  element 
cut  from  a  spherical  shell  by  two  adjacent  meridians  and  two  parallel  cir¬ 
cles.  The  notation  and  sign  convention  for  the  stress  resultants  conform 
to  that  of  Reference  11.  The  equations  of  equilibrium  of  the  element  are 
suinmarl/.ed  In  equation  (38)  anti  serve  to  determine  the  resultant  meridion¬ 
al  force  the  circumferential  force  A'#,  and  the  shearing  forces 

i\ "  A^.  All  of  the  forces  (tlat  is,  stress  resultants)  are  per  unit 
length  of  the  shell  section. 

3 

—  <A0  sin  4>)  +  Ne  cos  "  0; 

3  *Ne 

(A^  sin  <p)  +  -j£  +  cos  <p  -  0,  (38) 

£A0  +  Ag  +  F(4e)^  sin  <fi  -  0, 

where  P ^  is  the  pressure  distribution  normal  to  the  undeformed  shell 
surface.  A  solution  to  equation  (38)  is  sought  which  leads  to  finite 
stresses  at  the  shell  apex,  <p  =  0.  The  condition  that  the  resultants  of 
the  edge  tractions  on  the  base  circle  4>0  are  in  equilibrium  with  the  re¬ 
sultants  of  the  applied  forces  than  is  automatically  satisfied. 

Using  the  pressure  distribution  described  in  the  previous 
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section,  equation  (1)  ,  and  following  Reference  10,  a  solution  of  equation 
(38)  is  assumed  in  the  form 


and 


■  H  0  £  Su  (<£)  cos  n6, 

v  n»0  vn 

»(>  =  “  h  -  HP{<t>.6)> 


#4,6  “  Rq£o  st8n  W  sln  n 6' 


(39) 


Equation  (38)  is  then  reduced  to  a  pair  of  ordinary  differential  equations 
in  S^n  and  which  can  be  integrated  to'  yield 


4>n 


»  4  <0-  +  VJ 


S*8n  «  4  Wn  "  V0 


(40) 


Where,  letting 

C  ■  cos  0, 
S  a  sin  (/j, 
(/_ 


(1  +  C) 


ns  -  n  -  2 


U,«  W  -*x«  <°>l 


(41) 


and 


0  -  C)^'”'2  t*2n  <4>)-Kn], 


tin  W  f  <*>  ("  +  C>  <1  '  «"  S1""  <j& 


1  2n  W  f  nn  W  (n  -  0  (1  +  C)n  S 


n  ol-n 


d0, 


(48) 


J'l  2n  W  f°r  n  "  0,1 
arbitrary  for  n  >  2. 

V. 

The  integral  expressions,  equation  (42) ,  can  be  evaluated  once  the 
pressure  coefficients  Pn(<M  are  specified.  As  seen  from  equation  (2), 
these  will  generally  be  of  the  form 

t'n(4>)  =  2  j  +  2.  (sin  <£)"  t.  j  <  s. 

n  j* 0  n  J 

While  it  is  possible  to  evaluate  equation  (42)  for  the  j th  general  term, 
there  is  little  advantage  in  doing  so  for  small  A.  Explicit  results  for  . 
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N  =  3  are  given  below  for  the  Pn(<f>)  defined  in  equation  (2), 


— - - - 

Sj  0  K  ••  ' — ■  S 

2  4 


S#1  -  -  C(1  -  C)2  S"3 


J ,  .’13 

—  (2  +  0  +  —  (1  -  C)  (3C2  +  9C  i  8) 
3  5 


S^2  =  -  —  (2  +  B*>  +  A2  (1  +  C) 


-2 


S^3  - - (3  +  S2)  S  +  A3  S(l+C) 


-3 


=  0 


a2c 


-  —  -  A,  (1  +  cr2 


^  S  C -  A  3  (t+C’)-3, 

where,  as  before,  C  •*  cos  <p,  S  *  sin  0. 

We  observe  from  equation  (43)  or  equation  (42),  that  the  stress  func¬ 
tions  for  n  >  2  are  not  uniquely  defined.  This  results  from  the  fact  that 
the  pressure  terms  for  n  >  2,  equation  (1) ,  produce  no  resultant  force  or 
moment  on  the  shell,  and  it  would  be  necessary  to  specify  a  compatible 
distribution  of  stress  on  the  base  circle  in  order  to  establish  a  unique 
membrane  solution.  In  the  present  case,  of  course,  only  displacement  con¬ 
ditions  are  known  at  the  base  circle,  and  nothing  is  known  concerning  the 
stresses.  Thus,  the  membrane  solution  determined  is  unique  only  to  within 
an  arbitrary  system  of  self-equilibrated  stresses  applied  to  the  base  cir¬ 
cle.  The  way  out  is  to  consider  this  solution  as  the  first  step  in  a  two- 
step  procedure  leading  to  the  determination  of  the  bending  stresses  as  de¬ 
scribed  in  the  section  on  Strengths.  The  constants  of  integration,  An, 
appearing  in  equation  (43)  can  then  be  arbitrarily  determined  and  the  re¬ 
sulting  deformations  accounted  for  in  the  second  step  of  the  solution.  If 
one  considers  the  shell  to  be  fully  restrained  against  displacements  at 


§ 
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the  base  circle,  tho  Kn  are  most  conveniently  chosen  by  requiring  two  of 
the  three  components  of  the  base  displacements  to  be  zero,  thus  leaving 
only  one  displacement  to  be  eliminated  'by  base  moments  and  sheafs".  This 
evaluation  is  presented  in  the  following  section. 

SHELL  DISPLACEMENTS 

Following  Reference  12,  the  displacement  components  u,  u,  w  (measured, 
respectively,  in  the  positive  radial,  cucamferential  (d) ,  and  meridional 
(0)  directions)  can  be  related  to  the  stress  resultants  by  means  of  the 
following  system  of  partial  differential  equations. 


3b 

30  “ 

3w 

de 


W  cos  0  +  —  sin  4>  *  0  (N0  - N^)  sin  0, 
3n 

V  cos  0  -  sin  <(>  «  -  2/3  sin  0, 


(44) 


where, 


U  <=  —  (JV  ,  -  vNg)  +  ~  , 

1  +u  *  e  30  ’ 

„  ff(l  +  v) 

p’-rr~  • 


This  system  of  equations  can  be  integrated  in  a  fashion  similar  to  that 
employed  for  equation  (38).  Thus,  we  assume  a  solution  in  the  form 


V 

u  »  X  u  (<P)  cos  nff, 
»■  o  n 


x 

v  «  X  v  (0)  sin  rid, 

nB0  n 

a 

lo  ■  X  w  (0)  cos  riff. 

n=  o  n 


(48) 


Using  equations  (39)  and  (45)  and  after  considerable  manipulation,  the 
following  solutions  to  equation  (44)  are  obtained: 


--jjass1-"  t(i  -  a"  I/*  *-4,i  +  (1-0”  : 

«.  •  »-«•«*„-<>  («♦»>  S„  ♦  ».]  :  "  ’ 


U+^>  + 


where. 


Jxn  Kti  -  /  (SS0n  -  3S^fl  +  pn)  (1-0 
J2n  m  .  /  (2S*n  +  2  s^n  +  pn)  (1+C) 


n5  -  n  -1 


s  d. 


arbitrary  for  h  ■  0,1 


Jln  (0) 


for  ft  >  2 


■  arbitrary  for  tl  >  0, 

and  where  S^B  and  S^gfl  are  given  by  equations  (40) ,  (41)  and  (42)  and  pn 
is  determined  from  equation  (1). 

The  following  results  are  obtained  for  J  ln  and  J ln  for  H  *  3. 

a2c 

"  J2<>W  "  •  » 


<4 1  ^3 

J..(d>)  »  —  (2m“2+  2a-1  -  Im)  +  — (16m-2  + 

11  3  10 

+  16m-1  +  9m2  -  24m-  8lnm  +  45)  , 


^  i  ^3 

J,A<p)  K  —  Znm  +  —  (-9m2  +  12m  +  8 Inm  +  18), 
21  3  10 


JiiW  ‘  3  hm  +  —  . 


J.A4>)  •  -  K,  (1-20  m  -3—  (1-0  , 
3  ^  10 


‘'jjW  -  ~  A^m2 , 


JnW'T  » 
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where, 


ffi  ■  1  +  C  ■  1  +  cos  #« 

EVALUATION  QTThTEQRATIOM  CONSTANTS  ^ 

It  is  observed  from  the  results  of  the  preceding  two  sections  that  for 
terms  up  to  ?l  ■  A  (actually  A'  +  1  terms  since  n  •  0  is  counted),  there  is 
a  total  of  2A  +  2  integration  constants  consisting  of  A  -  1  constants  in  the 
stress  equation,  equation  (40),  and  N  +  3  constants  in  the  displacement 
equation,  (45) .  Since  there  are  3A’  +  3  components  of  displacement  to  be 
specified  at  the  base  circle,  it  is  evident  that  the  membrane  theory  can¬ 
not  possibly  serve  as  a  complete  solution  to  the  problem.  This,  of  course, 
was  stated  at  the  outset.  It  is  proposed  to  evaluate  these  constants  by 
setting  the  2N  «  2  components  of  the  l)  and  W  displacements  equal  to  zero 
at  the  base  circle,  $  ■  4>0>  This  will  provide  a  solution  for  a  shell  hav¬ 
ing  only  radial  displacements,  U  (<pQ)  «  u0,  at  the  base  circle.  The  com¬ 
plete  solution  for  the  fully  restrained  shell  then  can  be  obtained  by 
solving  the  general  bending  equations  for  a  system  of  edge  moments  and 
shears  which  remove  the  u0  displacements,  as  mentioned  previously. 

The  2A  +  2  equations  for  the  integration  constants  are  obtained  by 
setting  Vn(<pQ )  ■  W„(<p0)  *  O  in  equation  (46).  There  results 


for 

n 

n 

0,1 

O 

s: 

--4 

■ 

o 

for 

n 

> 

2, 

(40) 

\  “  J2n^ 

for 

n 

> 

o, 

where  the  J ln  contain  the  A  -  1  unknown  constants,  An(A  >  2).  Since 
vn  *  wn  "  0  on  4>  •  (pQ  ,  the  components  of  radial  displacement  at  the  base 
circle  are,  from  equation  (46)  , 
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[  (1  +  v)  S^n(tp0)  +  Pn(fp0)] 


(50) 
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For  A  ■ 
and  H.y 


3,  the  eight  integration  constants 
Letting  mQ  «  1  +  cos  these  are 


are  A  ^ ,  A  ^ .  L q . 
found  to  be, 
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A.  =  —  Mm.  +  —  (  -9  m„2  +  12  m„  +  8  In  mn+  15) , 

1  3  10  0  0  0 


Using  those  values,  the  stress  resultants  and  displacements  at  uny 
point  (<p  t  60  in  the  shell  can  bo  found  by  means  of  equations  (2),  (11)  , 
(311),  (13),  (45),  (40),  (48),  (50)  and  (51). 

HUMERI  CAL  RESULTS 

The  stress  resultants  and  displacements  for  A  »  3  were  programed  for 
computation  on  an  IISM-050  digital  computer.  Input  to  the  program  consists 
of  the  base  circle  angle  <£0 ,  Poisson's  ratio  v,  the  lead  parameters  40  ,  /!j, 
,12>  and  Ay,  and  information  as  to  the  number  and  location  of  points  in  the 
shell  at  which  stresses  and  displacements  are  to  be  determined.  Output 
consists  of  total  and  component  stresses,  deflections,  and  load  at  each 
point.  These  values  are  based  on  unit  stagnation  pressure,  i7,  and  unit 
shell  parameter  Et/lt2.  Actual  stresses  and  deflections  for  a  specific  load 
and  shell  are  obtained  by  a  simple  multiplication. 

Typical  resdlts  for  <£0  =  115°,  v  =  0.25,  ana  40  .  .  .  Ay  as  given  in 
Equation  (48)  are  shown  in  Figures  9  through  18.  Figures  9,  10,  and  11 
show  the  distribution  of  the  stress  resultants  on  a  series  of  parallel 
circles,  4>  ■  constant. 
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3  0  60  90  120  150  180 


0  (dtg) 

Figure  II.  Shear  Stress  Resultant  N4>6 

Inasmuch  as  the  membrane  stresses  are  not  an  adequate  representative 
of  total  stresses  in  the  vicinity  of  the  base  circle,  results  are  shown  only 
for  points  on  the  shell  lying  on  or  above  the  equatorial  plane,  that  is, 
for  <j>  >  90  degrees.  With  reference  to  Figures  a  and  10,  the  maximum  com¬ 
pressive  stresses  occur  at  the  stagnation  point,  the  circumferential  s tress 
being  three  and  a  half  times  as  great  as  the  meridional  stress.  The  maxi¬ 
mum  tensile  stresses  occur  in  the  vicinity  of  the  maximum  suction  pressure 
as  anticipated.  The  oscillatory  nature  of  the  stresses  in  the  back  por¬ 
tions  of  the  shell  probably  is  due  to  the  nature  of  the  approximate  pres¬ 
sure  distribution  (Figure  4).  The  maximum  combined  stresses  can  be  found 
in  a  straightforward  fashion  if  desired.  Figure  12  shows  the  stress  re¬ 
sultants  on  the  planes  A’  »  constant,  that  is,  on  sectors  of  the  shell 
whose  apex  Is  the  stagnation  point  and  whose  boundaries  coincide  with  the 
approximate  isobars  considered  (Figure  2)  .  The  contour  X  *=  0.91  defines 
a  circle  which  is  tangent  to  the  base  circle,  and,  as  Indicated  above, 
the  stresses  are  not  considered  to  be  meaningful  below  the  equatorial 
plane.  The  stresses  on  the  section  X  ■  0.99  are  believed  valid  at  all 
points;  the  situation  is  unclear  for  the  intermediate  section  X  »  0.95. 

The  components  of  displacement  are  shown  in  Figures  13,  14,  and  15. 

A  composition  of  displacements  showing  the  deformed  shape  of  the  shell  in 

>11 


the  planes  0  •  0  (vertical  plane  of  symmetry  containing  the  stagnation 
point)  $  ■  $o  “  liS  degrees  (the  base  plane)  ana  <p  ■  90  degrees  < 
equatorial  plane)  are-s-hown,'  respectively,-  HvFlgures- "TOJ1  IT  ana  iG ,  The 
displacement  pattern  of  Figure  17  suggests  the  nature  of  the  edge  moments 
and  shears  which  must  be  applied  to  effect  displacement  constraint  along 
the  base  circle. 


Figure  16.  Deformed  Shell  In  the  Plane,  Q  ■  0°,  1 80°(  u,  W  d  i  apl  acenenta) 
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load  parameters 
drag  coefficient 

lift  coefficient 

drag 

Young's  modulus  of  elasticity  of  shell  material 
integral  expressions,  Equations  (42)  and  (47) 
constants  of  integration 
lift 

resultant  circumferential  force  per  unit  length' ol‘ 
shell  cross  section 

resultant  meridional  force  per  unit  length  of  shell 
cross  section 

resultant  shearing  force  per  unit  length  of  shell 
cross  section 

pressure  distribution  parameter 

pressure  distribution  normul  to  the  undeformed 
shell  surface 

stagnation  pressure 
radius  <  r  shell 

meridional  and  shear  stress  components  defined  in 
Equation  (3'J) 

thickness  uf  shell 
radial  displacement  of  shell 
circumferential  displacement  of  shell 
meridional  uisplacement  of  shell 
coordinates,  see  Figure  3 

shell  parameter 

meridional  angle 
circumferential  angle 
angle  of  the  base  circle 
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